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EXHAUST GAS PURIFYING SYSTEM 
FOR INTERNAL COMBUSTION ENGINE 

BACKGROUND OF THE INVENTION 

[0001] The present invention relates to an exhaust gas 
5 purifying system of an internal combustion engine, and 
more particularly to a technique for preventing a 
particulate filter and a NOx trap catalyst from receiving 
excessive heat load during the recovery processing of 
these filter and catalyst, 

10 [0002] Particulate filters and NOx trap catalysts are 
common known as traps for removing specific contents from 
exhaust gas of an internal combustion engine. Each 
particulate filter has built in a filter element produced 
by molding ceramic into a honeycomb monolith. The filter 

15 element filters out particulates from exhaust gas. Each 
NOx trap catalyst changes its property according to the 
air/fuel ratio such as to remove NOx in exhaust gas by 
trapping NOx in the catalyst when the air/fuel ratio is 
lean. Such a NOx trap catalyst also traps sulfur content 

20 in exhaust gas in addition to NOx. These particulate 
filter and NOx trap catalyst are required to execute a 
recovery processing for recovering their performances 
when the accumulated quantity of eliminated objects such 
as particulates reaches a predetermine. If the engine is 

25 operated without executing the recovery processing of 

these filter and catalyst, there will cause an undesired 
increase of an engine back pressure and an undesired 
discharge of exhaust gas including NOx into atmosphere. 
Further, the NOx trap catalyst is required to execute 

30 recovery processing ( desulf urization recovery processing) 
for desulf urizing sulfur content trapped by NOx trap 
catalyst in addition to NOx. 
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[0003] Japanese Published Patent Application No- 2002- 
155793 discloses typical recovery processing of a 
particulate filter and a NOx trap catalyst wherein 
particulates trapped by the particulate filter are burnt 
5 by raising an exhaust gas temperature at a higher 

temperature than that during a normal operation, and NOx 
and sulfur content trapped by the NOx trap catalyst are 
discharged by temporally changing the air/fuel ratio. 
[0004] Japanese Published Patent Application No. 2000- 

10 179326 discloses a method of increasing an exhaust gas 
temperature by retarding a main injection timing, by 
executing a post injection, and by increasing a quantity 
of exhaust gas recirculation, for the recovery processing 
of a particulate filter and a NOx trap catalyst. 

15 SUMMARY OF THE INVENTION 

[0005] However, during the recovery processing of the 
particulate filter, the air /fuel ratio has been 
determined as a result of executing a post injection for 
reaching the exhaust gas temperature to a target 

20 temperature, and during the desulf urization recovery 

processing of the NOx trap catalyst, the air/fuel ratio 
has been determined as a result of supplying a reduction 
agent after raising the exhaust gas temperature. That is, 
no prior art has disclosed a technique of positively 

25 controlling an air/fuel ratio in the recovery processing. 
[0006] It is therefore an object of the present 
invention to provide an improved exhaust gas purifying 
system which is capable of recovering a particulate 
filter and a NOx trap catalyst without applying an 

30 excessive heat load to these filter and catalyst. 

[0007] An aspect of the present invention resides in 
an exhaust gas purifying system for an internal 
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combustion engine which comprises an exhaust gas 
purifying device which is disposed in an exhaust passage 
of the engine to remove specific content from exhaust gas 
and a control unit which is arranged to determine a 
5 recovery execution timing for executing recovery 

processing of recovering the exhaust gas purifying device 
from a specific content stacked state, to determine a 
target air/fuel ratio for executing the recovery 
processing, to determine a first engine controlled 

10 variable relating to an air/fuel ratio on the basis of 
the target air/fuel ratio, and to determine a second 
engine controlled variable relating to a combustion 
period at a value different from a value employed during 
normal processing when the recovery processing is 

15 executed - 

[0008] Another aspect of the present invention resides 
in an exhaust gas purifying system for an internal 
combustion engine, which comprises an exhaust gas 
purifying device disposed in an exhaust passage of the 

20 engine to remove specific content from exhaust gas and a 
control unit which is arranged to determine whether 
recovery processing for recovering the exhaust gas 
purifying device as to accumulated specific contents in 
the exhaust gas purifying device is executed, and to 

25 increase an exhaust gas temperature at a temperature 
higher than an exhaust gas temperature during normal 
processing, by setting an air/fuel ratio at a target 
air/fuel ratio and by controlling the a combustion period 
while maintaining the air/fuel ratio at the target 

30 air/fuel ratio when the recovery processing is executed, 
[0009] A further aspect of the present resides in a 
method of executing recovery processing of an exhaust gas 
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purifying disposed in an exhaust passage of an internal 
combustion engine. The method comprises an operation of 
determining a recovery execution timing for recovery 
processing of recovering the exhaust gas purifying device 
5 from a specific content stacked state, an operation of 
setting a target air/fuel ratio for executing the 
recovery processing, an operation of setting a first 
engine controlled variable relating to an air/fuel ratio 
on the basis of the target air/fuel ratio, and an 
10 operation of setting a second engine controlled variable 
relating to a combustion period at a value different from 
a value employed during normal processing when the 
recovery processing is executed . 

[0010] The other objects and features of this 
15 invention will become understood from the following 

description with reference to the accompanying drawings. 
BRIEF DESCRIPTION OF THE DRAWINGS 

[0011] Fig. 1 is a view showing a direct injection 
type diesel engine provided with an embodiment of an 
20 exhaust gas purifying system according to the present 
invention . 

[0012] Fig. 2 is a block diagram showing an electronic 
control unit of the embodiment according to the present 
invention. 

25 [0013] Fig. 3 is a flowchart of a mode decision value 
setting routine. 

[0014] Fig. 4 is a flowchart of a target acceleration 
request injection quantity calculation routine. 
[0015] Fig. 5 is a map for obtaining an target 
30 acceleration request injection quantity. 

[0016] Fig. 6 is a flowchart of an intake system 
response time constant calculation routine. 
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[0017] Fig. 7 is a map for obtaining a volumetric 
efficiency basic value, 

[0018] Fig. 8 is a table for obtaining a volumetric 
efficiency correction value. 
5 [0019] Fig. 9 is a flowchart of a cylinder intake air 
quantity calculation routine. 

[0020] Fig. 10 is a conversion table between a voltage 
and an intake air quantity. 

[0021] Fig. 11 is a flowchart of an exhaust gas flow 
10 rate calculation routine. 

[0022] Fig. 12 is a flowchart of an EGR rate 
calculation routine . 

[0023] Fig. 13 is a flowchart of a turbine nozzle 
opening calculation routine 
15 [0024] Fig. 14 is a flowchart of an EGR gas flow 
velocity calculation routine. 

[0025] Fig. 15 is a map for obtaining an EGR gas flow 

velocity basic value. 

[0026] Fig. 16 is a map for obtaining an EGR gas flow 
20 velocity correction value. 

[0027] Fig. 17 is a flowchart of a recovery mode 

target excess air ratio calculation routine. 

[0028] Fig. 18 is a map for obtaining a target excess 

air ratio basic value. 
25 [0029] Fig. 19 is a flowchart of an excess air ratio 

calculation routine. 

[0030] Fig. 20 is a conversion table between a pump 
current and the excess air ratio. 

[0031] Fig. 21 is a flowchart of a torque correction 
30 coefficient calculation routine. 

[0032] Figs. 22A and 22B are maps for obtaining torque 
correction coefficients . 
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[0033] Fig. 23 is a flowchart of a target intake air 
quantity calculation routine. 

[0034] Fig. 24 is a flowchart of a target fuel 
injection quantity calculation routine. 
5 [0035] Fig. 25 is a flowchart of an intake throttle 
valve opening calculation routine. 

[0036] Fig. 26 is a table for obtaining a maximum 
working gas quantity. 

[0037] Fig. 27 is a table for obtaining an intake air 
10 quantity ratio. 

[0038] Fig. 28 is a conversion table between an 
opening area and a valve opening . 

[0039] Fig. 29 is a flowchart of a target EGR rate 
basic value calculation routine. 
15 [0040] Fig. 30 is a map for obtaining a target EGR 
rate basic value. 

[0041] Fig. 31 is a flowchart of a target EGR rate 
calculation routine. 

[0042] Fig. 32 is a flowchart of a target EGR gas 
20 quantity calculation routine. 

[0043] Fig. 33 is a flowchart of a target EGR valve 
opening calculation routine. 

[0044] Fig. 34 is a conversion table between a valve 
opening and number of steps. 
25 [0045] Fig. 35 is a flowchart of a target turbine 
nozzle opening calculation routine. 

[0046] Fig. 36 is a map for obtaining a turbine nozzle 
opening basic value. 

[0047] Fig. 37 is a map for obtaining a turbine nozzle 
30 opening correction value. 

[0048] Fig. 38 is a flowchart of a target turbine 
nozzle opening delay compensation routine. 
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[0049] Fig. 39 is a table for obtaining an exhaust 
system response time constant. 

[0050] Fig. 40 is a table for obtaining an advance 
compensation coefficient . 
5 [0051] Fig. 41 is a flowchart of a target duty ratio 
calculation routine. 

[0052] Fig. 42 is a conversion table between a nozzle 
opening and a duty ratio. 

[0053] Fig. 43 is a flowchart of a main injection 
10 timing calculation routine. 

[0054] Fig. 44 is a map for obtaining a main injection 
timing basic value. 

[0055] Fig. 45 is a flowchart of an air/fuel ratio 
feedback control routine. 
15 [0056] Figs. 46A # 46B and 46C are tables for obtaining 
compensation gains . 

[0057] Fig. 47 is a flowchart of a target exhaust gas 
temperature calculation routine. 

[0058] Fig. 48 is a flowchart of an exhaust gas 
20 temperature calculation routine. 

[0059] Fig. 49 is a conversion table between a voltage 

and a temperature as to the exhaust gas temperature. 

[0060] Fig. 50 is a flowchart of an exhaust gas 

temperature feedback control routine. 
25 [0061] Figs. 51A, 51B and 51C are tables for obtaining 

compensation gains . 

[0062] Figs. 52A, 52B and 52C are graphs showing a 
relationship of the exhaust gas temperature, a CO 
discharge quantity and a HC discharge quantity relative 
30 to the air/fuel ratio, a flowchart of a target exhaust 
gas temperature calculation routine. 
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[0063] Fig. 53 is a graph showing a relationship 
between a heated temperature of a NOx trap catalyst and a 
NOx conversion ratio. 

[0064] Figs. 54A and 54B is graphs showing 
5 relationships of the exhaust gas temperature and a 

particulate combustion speed relative to the air/fuel 
ratio . 

[0065] Figs. 55A, 55B and 55C are views for explaining 
a malfunction of a diesel particulate filter. 

10 DETAILED DESCRIPTION OF THE INVENTION 

[0066] Referring to the drawings, there is discussed 
an embodiment of an exhaust gas purifying system 
according to the present invention. Fig. 1 shows a 
direct injection diesel engine 1 which employs the 

15 exhaust gas purifying system according to the present 
invention . 

[0067] An air cleaner 12 for removing dust particles 
from intake fresh air is attached to an inlet of an 
intake passage 11. An airflow meter 13 is disposed 

20 downstream of air cleaner 12 and measures an airflow rate. 
The air passed through air cleaner 12 and airflow meter 
13 is flowed into a collector 14 and is distributed into 
cylinders through an intake manifold. 
[0068] A nozzle-variable type turbocharger 15 is 

25 attached to engine 1, and more specifically a compressor 
section 15a of turbocharger 15 is disposed upstream of 
collector 14. An intercooler 16 is disposed between 
compressor section 15 and collector 14 to cool the intake 
air compressed by turbocharger 14. An intake throttle 

30 valve 14, through which intake air flow rate is 

controlled, is provided upstream of collector 14. A 
swirl control valve 18 for controlling gas flow in each 
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cylinder is provided at each intake port for each 
cylinder. An electronic control unit (ECU) 61 outputs 
control signals to intake throttle valve 17 and swirl 
control valve 18 , respectively. 
5 [0069] Fuel injectors 21 for the respective cylinders 
are fixed to a cylinder head of an engine body so that an 
injecting portion of each fuel injector 21 is faced with 
a combustion chamber upper portion of each cylinder. A 
fuel system of engine 1 comprises a common rail 22 so 

10 that fuel fed by a fuel pump 23 is controlled at a 

predetermined pressure and is supplied through common 
rail 22 to each fuel injector 21. Each fuel injector 21 
executes fuel injection in response to the signal from 
ECU 61. The fuel injection by each fuel injector 21 is 

15 constructed by several time injections. Each fuel 

injector 21 executes a pilot injection at a moment before 
a main injection is executed in addition to the main 
injection. This pilot injection suppresses the 
generation of particulates and reduces the level of 

20 combustion noise. 

[0070] A NOx trap catalyst 32, which traps NOx or 
reduces and desorbs the trapped NOx according to the 
air/fuel ratio of exhaust gas, is disposed downstream of 
an exhaust manifold of engine 1 . A diesel particulate 

25 filter 33 functioning as a particulate filter is disposed 
downstream of NOx trap catalyst 32. Under a normal lean 
operating condition, NOx and particulates in exhaust gas 
are removed from the exhaust gas by NOx trap catalyst 32 
and diesel particulate filter 33. Under the normal lean 

30 operating condition, NOx trap catalyst 32 traps sulfur in 
the exhaust gas in addition to NOx. 
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[0071] An EGR conduit 41 connects exhaust passage 31 
and intake passage 11. An EGR valve 42 is disposed in 
EGR conduit 41. By controlling an opening of EGR valve 
42 according to the control signal of ECU 61, a proper 
5 quantity of exhaust gas according to the opening degree 
of EGR valve 42 is recirculated (returned) to intake 
passage 11. An EGR gas cooler 43 is disposed upstream of 
EGR valve 4 2 to cool EGR gas. 

[0072] A turbine 15b of turbocharger 15 is disposed 

10 between a portion connected to EGR conduit 42 and NOx 

trap catalyst 32 in exhaust passage 31. A nozzle opening 
of turbine 15b is controlled by an actuator 51 which 
moves a variable vane of turbine 15b according to a 
signal outputted from ECU 61. 

15 [0073] The exhaust gas purifying system comprises NOx 
trap catalyst 32, diesel particulate filter 33, ECU 61 
having a recovery function of catalyst 32 and filter 33, 
and sensors. These sensors includes airflow meter 13, a 
sensor 71 for detecting a temperature Tw of engine 

20 coolant, a sensor 72 for detecting an excess air ratio 
lamb of exhaust gas upstream of NOx trap catalyst 32, a 
sensor 73 for detecting an exhaust gas temperature Texh 
of exhaust gas upstream of diesel particulate filter 33, 
and a sensor 74 for detecting a pressure difference APdpf 

25 between pressures at upstream side and downstream side of 
diesel particulate filter 33. 

[0074] Fig. 2 is a block diagram showing functions of 
ECU 61. In Fig. 2, a module Ml performs a function of 
setting a mode decision value ATSstate to changeover an 
30 operation mode between a normal mode and a recovery mode. 
The recovery mode includes a desulf urization recovery 
mode and a filter recovery mode. A module M2 performs a 
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function of calculating an inner state quantity of engine 
1, such as a cylinder intake air quantity Qac, and an EGR 
rate Regr. A module M3 performs a function of 
calculating a target excess air ratio tlamb according to 
5 mode decision value ATSstate. A module M4 performs a 

function of calculating an actual excess are ratio lamb. 
A module M5 performs a function of calculating a target 
EGR rate, a target turbine opening, a target intake 
throttle opening and a target fuel injection quantity to 

10 achieve target excess air ratio tlamb. A module M6 

performs functions of calculating a difference between 
target excess air ratio tlamb and actual excess air ratio 
lamb and of calculating a feedback correction quantity 
for one of fuel injection quantity and intake throttle 

15 valve opening on the basis of the obtained difference so 
as to bring actual excess air ratio lamb closer to target 
excess air ratio tlamb. A module M7 performs a function 
of correcting a main injection timing so as to achieve a 
target exhaust gas temperature tTexh according to mode 

20 decision value ATSstate when the recovery mode is 

selected. Further, module M7 performs a function of 
correcting a pilot injection timing to suppress the 
generation of smoke and to reduce the generation of 
combustion noise. The main injection timing and the 

25 pilot injection timing correspond to a second engine 

controlled variable. The pilot injection quantity may be 
employed as second engine controlled variable. A module 
M8 performs a function of calculating exhaust gas 
temperature Texh on the basis of the signal outputted 

30 from sensor 73. A module M9 performs a function of 

correcting the main injection timing so as to decrease a 
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difference between target exhaust gas temperature tTexh 
and exhaust gas temperature Texh. 

[0075] Hereinafter, there are discussed operations of 
ECU 61 corresponding to modules Ml through M9 . 
5 [0076] First there is discussed a manner of setting of 
mode decision value ATSstate. Fig. 3 is a flowchart of a 
mode decision value setting routine and is started at 
module Ml of ECU 61. 

[0077] At step S101 ECU 61 reads coolant temperature 
10 Tw, exhaust gas flow rate Qexh, an engine speed Ne # 
exhaust gas temperature Texh. At step S102 ECU 61 
determines whether or not coolant temperature Tw is 
higher than or equal to a predetermined temperature Twl . 
When the determination at step SI is negative, that is, 
15 coolant temperature Tw is lower than predetermined 
temperature Twl, the routine proceeds to step S103 
wherein ECU 1 sets mode decision value ATSstate at 0 
( ATSstate=0 ) . Thereafter, the routine returns to a start 
block. When the determination at step SI is affirmative, 
20 that is, coolant temperature Tw is higher than or equal 

to predetermined temperature Twl, the routine proceeds to 
step S104 wherein ECU 61 sets mode decision value 
ATSstate at 1 ( ATSstate=l ) . 

[0078] At step S105 subsequent to the execution of 
25 step S103, ECU 61 calculates a NOx exhaust quantity NOX 
per unit time in exhaust gas, on the basis of exhaust gas 
flow rate Qexh. At step S106 ECU 61 obtains an integral 
of NOx exhaust quantity NOX and stores the obtained 
integral as a NOx trap quantity 2 NOX, which is stored in 
30 NOx trap catalyst 32, in a memory of ECU 61. At step 
S107 ECU 61 calculates an integrating value of engine 
speed Ne and stores the obtained value as a sulfur trap 
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quantity ESOX, which is stored in NOx trap catalyst 32, 
in the memory of ECU 61. Further, at step S107 ECU 61 
stores the obtained value as a particulate accumulated 
quantity EPM, which is stored in diesel particulate 
5 filter 33, in the memory of ECU 61. 

[0079] At step S108 ECU61 determines whether or not 
NOx trap quantity ENOX is greater than or equal to a 
predetermined quantity ENOX1. When the determination at 
step S108 is affirmative ( ENOX^ ENOX1 ) , the program 
10 proceeds to step S109 wherein ECU 61 sets mode decision 

value ATSstate at 2 ( ATSstate=2 ) . When the determination 
at step S108 is negative ( ENOX <EN0X1), the program 
proceeds to step S110. 

[0080] At step S110 ECU 61 determines whether or not 
15 exhaust gas temperature Texh is higher than or equal to a 
predetermined temperature Texhl . When the determination 
at step S110 is affirmative ( Texh^Texhl ) , the program 
proceeds to step Sill. When the determination at step 
S110 is negative ( Texh<Texhl ) , the program proceeds to a 
20 return block to return the present routine. 

[0081] At step Sill ECU 61 determines whether or not 
sulfur trap quantity ESOX is greater than or equal to a 
predetermined quantity ESOX1. When the determination at 
step Sill is affirmative ( £ SOX^ £ SOX1 ) , the program 
25 proceeds to step S112 wherein ECU 61 set mode decision 

value ATSstate at 3 ( ATSstate=3 ) . When the determination 
at step Sill is negative ( £ SOX<£ SOX1 ) , the program 
proceeds to step S113. 

[0082] At step S113 ECU 61 determines whether or not 
30 particulate accumulated quantity EPM is greater than or 



P04NM- 007/02 -00560 

-14- 

equal to a predetermined quantity 2PM1. When the 
determination at step S113 is affirmative ( 2 PM^ E PM1 ) , 
the program proceeds to step S114 wherein ECU 61 sets 
mode decision value ATSstate at 1 ( ATSstate=l ) . 
5 Thereafter the present routine is returned. When the 
determination at step S113 is negative (£PM<2PM1), the 
program proceeds to the return block to return the 
present routine, 

[0083] When mode decision value ATSstate is set at 3 

10 or 4, ECU 61 selects an exhaust gas temperature rising 
mode and executes a control for rising the exhaust gas 
temperature under a condition that the excess air ratio 
is set at a stoichiometric air/fuel ratio or neighborhood 
thereof. After the exhaust gas temperature reaches the 

15 target temperature, when mode decision value ATSstate is 
set at 3 (ATSstate=3) , ECU 61 selects a desulf urization 
mode and executes a control for discharging sulfur 
trapped in NOx trap catalyst 32 by varying the excess air 
ratio to a rich side. When mode decision value ATSstate 

20 is set at 4 ( ATSstate=4 ) , ECU 61 selects a filter 
recovery mode and executes a control for burning 
particulates accumulated in diesel particulate filter 33 
by varying the excess air ratio to a lean side. When 
exhaust gas temperature Texh reaches a second 

25 predetermined temperature Texh2 higher than predetermined 
temperature Texhl during the recovery processing, ECU 61 
selects a malfunction avoidance mode, lowers exhaust gas 
temperature Texh by varying the excess air ratio to the 
lean side, and suspends the recovery processing, in order 

30 to prevent the functional deterioration of NOx trap 

catalyst 32 or burnout of diesel particulate filter 33. 
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[0084] Hereinafter, there is discussed calculations of 
inner state quantities. Fig. 4 is a flowchart for a 
calculation routine of a target accelerating demand fuel 
injection quantity, which is executed by module M2 of ECU 
5 61. 

[0085] At step S201 ECU 61 reads engine speed Ne and 
control lever opening APO. At step S202 ECU 61 retrieves 
(calculates) an acceleration demand fuel injection basic 
value Mqdrv from a map shown in Fig . 5 , engine speed Ne 

10 and control lever opening APO. At step S203 ECU 61 

calculates an idling speed correction quantity Qfisc. At 
step S204 ECU 61 calculates a target acceleration demand 
fuel injection quantity Qfdrv by adding the obtained 
speed correction quantity Qfisc to acceleration demand 

15 fuel injection basic value Mqdrv ( Qf drv=Mqdrv+Qf isc ) . 

[0086] Fig. 6 is a flowchart for a calculation routine 
of an air intake system response time constant, which is 
executed by module M2 of ECU 61. 

[0087] At step S211 ECU 61 reads engine speed Ne, 
20 target acceleration demand fuel injection quantity Qfdrv, 
an intake manifold pressure Pint and a target EGR rate 
Megrdn-i. Herein, a first-order delayed value Megrd of 
target EGR rate Megr is approximated as an actual EGR 
rate, and reference n-1 denotes that the value with this 
25 reference n-1 was obtained in the previous routine. 

[0088] At step S212 ECU 61 retrieves (calculates) a 
volumetric efficiency basic value Kinb from a map shown 
in Fig. 7 on the basis of engine speed Ne and target 
acceleration demand fuel injection quantity Qfdrv, and 
30 further retrieves (calculates) a volumetric efficiency 
correction value Kinh from a map shown in Fig. 8 on the 
basis of intake manifold pressure Pint. 
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[0089] AT step S213 ECU 61 calculates a volumetric 
coefficient Kin on the basis of volumetric efficiency 
basic value Kinb, volumetric efficiency correction value 
Kinh and target EGR rate Megrd n _i from the following 
5 expression (1). 

Kin= Kinb X Kinh X (i/(i+Mergd n .i) ) ---(1) 
[0090] At step S214 ECU 61 calculates an air intake 
system response time constant Kkin by multiplying 
volumetric coefficient Kin by a volumetric ratio KVOL# 

10 (Kkin=KinXKVOL#) . Volumetric ratio KVOL# is a ratio of a 
stroke volume Vc of pistons Vc and a volume Vm of the 
intake manifold including the collector 14 (KVOL#=Vc/Vm) . 
[0091] Fig, 9 is a flowchart for a calculation routine 
of an cylinder intake air quantity, which is executed by 

15 module M2 of ECU 61. 

[0092] At step S221 ECU 61 reads output AFM of airflow 
meter 13 , engine speed Ne, air intake system response 
time constant Kkin. At step S222 ECU 61 obtains an 
intake air quantity Qas by converting airflow meter 

20 output AFM using a table shown in Fig. 10. At step S223 
ECU 61 obtains a basic value QasO by executing a weighted 
average process of intake air quantity Qas. At step S224 
ECU 61 obtains a per-cylinder per-stroke intake air 
quantity QacO from the following expression (2). 

25 Qac0= (AasO/Ne) XKCON# (2) 

where KCON# is a unit conversion coefficient. 

[0093] At step S225 ECU 61 calculates a collector 
inlet intake air quantity Qacn by executing n- times delay 
processing of per-cylinder per-stroke intake air quantity 

30 QacO (Qacn=Qac0 n -k) • At step S226 ECU 61 calculates a 
cylinder intake fresh air quantity Qac by executing a 
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delay processing of collector inlet intake air quantity 
Qacn using the following expression (3). 

Qac= Qac n -i X ( 1 - Kkin ) +Qacn X Kkin (3) 

[0094] Fig. 11 is a flowchart for a calculation 
5 routine of an exhaust gas flow rate , which is executed 
by module M2 of ECU 61. 

[0095] At step S231 ECU 61 reads cylinder intake air 
quantity Qac, EGR gas quantity (Qec=tQec0), target 
accelerating demand fuel injection quantity Qfdrv and 

io engine speed Ne. At step S232 ECU 61 obtains a unit time 
quantity Qf of target accelerating demand fuel injection 
quantity Qfdrv using the following expressions (4). 

Qf =Qf drvX Ne/KCON# ( 4 ) 

[0096] At step S233 ECU 61 obtains a unit time 

15 quantity Qa of cylinder intake air quantity Qac using the 
following expression (5). 

Qa=Qac X Ne/KCON# ( 5 ) 

[0097] At step S234 ECU 61 obtains a unit time 
quantity Qe of EGR gas quantity Qec using the following 

20 expression (6). 

Qe=QecXNe/KCON# (6) 

[0098] At step S235 ECU 61 calculates an exhaust gas 
flow rate on the basis of the obtained quantities Qf # Qa 
and Qe using the following expression (7). 

25 Qexh=Qa+Qe+Qf XGKQF# (7) 

[0099] Fig. 12 is a flowchart for an EGR rate 
calculation routine, which is executed by module M2 of 
ECU 61. 

[00100] As discussed above, target EGR rate Megr is 
30 approximated by a first -order delay value of actual EGR 
rate. Accordingly at step S241 ECU 61 reads target EGR 
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rate Megr n . lf target EGR gas quantity tQecd n -i and cylinder 
intake air quantity Qac. 

[00101] At step S242 ECU 61 obtains the first-order 
delay value Megrd by executing a first-order delay 
5 processing of target EGR rate Megr n -i using the following 
expression (8), and stores the obtained first-order delay 
value Mefrd. 

Megrd=(l-TCECR#) XMegrd n _i+TCEGR# XMegr n _! 

---(8) 

10 [00102] At step S243 ECU calculates an EGR rate Regr by 
dividing target EGR gas quantity tQecd n _i by cylinder 
intake air quantity Qac as shown by the following 
expression ( 9 ) . 

Regr=tQece n _i/Qac (9) 

15 [0103] Fig. 13 is a flowchart for a calculation 

routine of a turbine nozzle opening, which is executed by 
module M2 of ECU 61. 

[0104] Target turbine nozzle opening Trav is 
approximated by a first -order delay value of an actual 

20 turbine nozzle opening. Accordingly at step S251 ECU 61 
reads a target turbine nozzle opening Travff n -i. At step 
S252 ECU 61 obtains turbine nozzle opening Rvgt by 
executing a first -order delay processing of target 
turbine nozzle opening Travff n -i using the following 

25 expression (10), and stores the obtained turbine nozzle 
opening Rvgt . 

Rvgt=(l-TCVGT#) XRvgtn-i+TCVGTtfXTravffn-i 

---(10) 

[0105] Fig. 14 is a flowchart for a calculation 
30 routine of an EGR gas flow velocity, which is executed by 
module M2 of ECU 61. 
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[0106] An EGR gas flow velocity Cqe is obtained on the 
basis of intake manifold pressure Pint, an exhaust 
manifold pressure Pexh and an exhaust gas gravity, using 
the following expression (11). 

5 Cqe=V(2px(Pexh-Pint) ) (11) 

However, it is difficult to accurately measure intake 
manifold pressure Pint and exhaust manifold pressure Pexh. 
Accordingly, EGR gas flow velocity Cqe is estimated by 
the following method. 

10 [0107] At step S261 ECU 61 reads EGR gas quantity Qec 
(=tQecd), intake air quantity Qacn, turbine nozzle 
opening Rvgt and intake throttle valve opening TVO. At 
step S262 ECU 61 retrieves (calculates) a flow velocity 
basic value CqeO from a map shown in Fig. 15 using gas 

15 quantity Qec and intake throttle valve opening TVO. At 
step S263 ECU 61 retrieves a flow velocity correction 
value Kcqe from a map shown in Fig. 16 using intake air 
quantity Qacn and turbine nozzle opening Rvgt. At step 
S264 ECU 61 calculates EGR gas flow velocity Cqe by 

20 multiplying flow velocity basic value CqeO by flow 
velocity correction value Kcqe (Cqe=Cqe0 XKcqe) . 
[0108] Herein, there is discussed a setting of a 
target air /fuel ratio. Fig. 17 is a block diagram 
showing a calculation routine of a target excess air 

25 ratio, which is executed by module M3 of ECU 61. 

[0109] At steps S301 and S302 ECU 61 reads mode 
decision value ATSstate and selects a map corresponding 
to mode decision value ATSstate. Further, ECU 61 
retrieves a target excess air ratio basic value TlambO 

30 according to the operation mode from the selected map. 
[0109] More specifically, when ATSstate=0, ECU 61 
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searches a low temperature target A map and sets target 
excess air ratio basic value TlambO at 1 indicative of a 
stoichiometric air/fuel ratio. When ATSstate=l, ECU 61 
searches a normal target A map shown in Fig. 18 and sets 
5 target excess air ratio basic value TlambO at 1.4 or more 
indicative of a lean state. When ATSstate=2, ECU 61 sets 
target excess air ratio basic value TlambO at 0.9 
indicative of a rich state. When ATSstate=3, ECU 61 
searches a desulf urization mode target A map and sets 

10 target excess air ratio basic value TlambO at 0.99 

indicative of the rich state. When ATSstate=4, ECU 61 
searches a filter recovery mode target A map and sets 
target excess air ratio basic value TlambO at 1.2 
indicative of the lean state. 

15 [0110] When ATSstate=3 or 4 , an exhaust gas rising 

mode is executed before the desulf urization mode or the 
filter recovery mode are executed. During this exhaust 
gas rising mode, ECU 61 sets target excess air ratio 
basic value TlambO at 1 indicative of a stoichiometric 

20 air/fuel ratio. When the processing of the 

desulf urization or filter cleaning is suspended due to 
the excessive rising of the exhaust gas temperature, ECU 
61 sets target excess air ratio basic value TlambO at 1.3 
or more. That is. During the malfunction avoidance mode, 

25 ECU 61 sets target excess air ratio basic value TlambO at 
a value greater than that during the filter recovery mode. 
[0111] At step S303 ECU 61 executes a delay processing 
of target excess air ratio basic value TlambO using the 
following expression (12) employing intake system 

30 response time constant Kkin and obtains a target excess 
air ratio Tlamb. 
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Tlamb=Tlamb n -i X (1-Kkin) + TlambO X Kkin (12) 

[0112] There is discussed a calculation of the excess 
air ratio. Fig. 19 is a flowchart for an excess air 
ratio calculation routine, which is executed by module M4 
5 of ECU 61. 

[0113] At step S401 ECU 61 reads a pump current ip 
from sensor 72. At step S402 ECU 61 obtains excess air 
ratio lambO from a table shown in Fig. 20 using pump 
current ip. At step S403 ECU 61 executes a weighted 

10 average processing of excess air ratio lambO and sets the 
obtained value as excess air ratio lamb. 
[0114] There is discussed a setting of an engine 
controlled variable. Fig. 21 is a flowchart for a 
calculation routine of a torque correction coefficient, 

15 which is executed by module M5 of ECU 61. ECU 61 

determines a torque correction coefficient Ka according 
to target excess air ratio Tlamb and main injection 
timing MITf and uses the obtained torque correction 
coefficient Ka in a target intake air quantity 

20 calculation routine and a target fuel injection 
calculation routine. 

[0115] At step S501 ECU 61 reads target excess air 
ratio Tlamb, engine speed Ne and main injection quantity 
MITf. At step S502 ECU 61 retrieves a first torque 

25 correction coefficient KaLAMB from a map shown in Fig. 
22A with reference to target excess air ratio Tlamb and 
engine speed Ne and retrieves a second torque correction 
coefficient KaMIT from a map shown in Fig. 22B with 
reference to main fuel injection timing MITf and engine 

30 speed Ne. The first torque correction coefficient KaLAMB 
is set to adapt to a change of target excess air ratio 
Tlamb during the recovery mode and is set at a value. 
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which is greater than 1 and increases as target excess 
air ratio Tlamb is decreased, when target excess air 
ratio Tlamb is smaller than 1.4. Further, the first 
torque correction coefficient KaLAMB is set at 1 when 
target excess air ratio Tlamb is greater than or equal to 
1.4. On the other hand, the second torque correction 
coefficient KaMIT is set to adapt to a change of main 
injection timing MITf during the recovery mode, and is 
set at a value, which is greater 1 when main injection 
timing MITf is retarded relative to a normal timing MITO 
and which increases as the degree of the retard of main 
injection timing MITf increases. Second torque 
correction coefficient KaMIT is normally set at 1 . 
[0116] At step S503 ECU 61 obtains torque correction 
coefficient Ka by multiplying first torque correction 
coefficient KaLAMB and second torque correction 
coefficient KaMIT (Ka=KaLAMB X KaMIT). 

[0117] Fig. 23 is a flowchart for a target intake air 
quantity calculation routine, which is executed by module 
M5 of ECU 61. 

[0118] At step S511 ECU 61 reads target excess air 
ratio Tlamb, target acceleration demand injection 
quantity Qfdrv and torque correction coefficient Ka. At 
step S512 ECU 61 calculates a target intake air quantity 
basic value tQacO from the following expression (13) on 
the basis of target excess air ratio Tlamb, target 
acceleration demand injection quantity Qfdrv and torque 
correction coefficient Ka. 

tQacO=Tlamb X Qfdrv X Blamb# X Ka ---(13) 
where Blamb3 is a stoichiometric air/fuel ratio 
corresponding value (14.7). 
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[0119] At step S513 ECU 61 executes a weighted average 
processing of target intake air quantity basic value 
tQacO and sets the obtained value as target intake air 
quantity tQac . 
5 [0120] Fig- 24 is a flowchart for a target fuel 

injection quantity calculation routine, which is executed 
by module M5 of ECU 61. 

[0121] At step S521 ECU 61 reads target excess air 
ratio Tlamb, intake air quantity Qac, target accelerating 

10 request injection quantity Qfdrv, torque correction 

coefficient ka and mode decision value ATSstate. At step 
S522 ECU 61 determines whether or not mode decision value 
ATSstate is one of 0, 2 and 3. When the determination at 
step S522 is affirmative, that is, when mode decision 

15 value ATSstate is one of 0, 2 and 3, the air/fuel ratio 
is controlled at a rich state or stoichiometric state, 
and therefore engine torque is mainly dependent on intake 
fresh air. Accordingly the program proceeds to step S523 
wherein ECU 61 calculates target fuel injection quantity 

20 tQf using the following expression (14) on the basis of 
intake air quantity Qac. 

tQf =Qac/( Tlamb X Blamb# ) X Ka (14) 

On the other hand, when the determination at step S522 is 
negative, that is, when mode decision value ATSstate is 

25 neither of 0, 2 nor 3, the air/fuel ratio is controlled 
at lean state, and therefore the engine torque is mainly 
determined by the fuel injection quantity. Accordingly, 
the program proceeds to step S524 wherein ECU 61 
calculates target fuel injection quantity tQf using the 

30 following expression (15) on the basis of target 
accelerating request injection quantity Qfdrv. 

tQf=Qfdrv X Ka (15) 
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[0122] Fig. 25 is a flowchart for an intake throttle 
valve opening calculation routine, which is executed by 
module M5 of ECU 61. 

[0123] At step S531 ECU 61 reads engine speed Ne, 
5 target EGR rate Megr and target intake air quantity tQac. 
At step S523 ECU 61 retrieves a maximum working gas 
quantity Qgmax from a table shown in Fig. 2 6 with 
reference engine speed Ne. At step S533 ECU 61 
calculates a target working gas quantity ratio tQhO from 
10 the following expression (16) on the basis of target 
intake air quantity tQac. 

tQhO=tQac X (1+Megr)/VCE# /Qgmax (16) 

where VCE# is a stroke volume of piston. 

[0124] At step S534 ECU 61 obtains a target air flow 

15 rate tDNV through a conversion of target working gas 

quantity ratio tQhO using a table shown in Fig. 27. At 
step S535 ECU 61 calculates a target opening area basing 
value tAtvob from the following expression (17) on the 
basis of target air flow rate tDNV and engine speed Ne. 

20 tAtvob=tDNV X Ne X VOL# (17) 

[0125] At step S536 ECU 61 calculates a target intake 
throttle value opening area tAtvo from the following 
expression (18) on the basis of target opening area 
basing value tAtvob and target EGR rate Megr. 

25 tAtvo=tAtvob X 1/(1+Megr) (18) 

where tAtvo is a value obtained by correcting target 
opening area basic value tAtvob, which is a target 
opening area with respect to the total working gas , by 
target EGR rate Megr. At step S537 ECU 61 obtains intake 

30 throttle valve opening ETC through a conversion of target 
intake throttle valve opening area tAtvo using a table 
shown in Fig. 28. 
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[0126] Fig. 29 is a flowchart for a calculation 
routine of a target EGR rate basic value, which is 
executed by module M5 of ECU 61. 

[0127] At steps S541 and S542 ECU 61 reads mode 
5 decision value ATSstate and selects a map corresponding 
to mode decision value ATSstate. Further, ECU 61 
retrieves target EGR rate basis value MegrO according to 
the operation mode from the selected map. 
[0128] More specifically, when ATSstate=l, ECU 61 

10 searches a standard map shown in Fig. 30 and sets normal 
value as target EGR rate basis value MegrO. When 
ATSstate=0, ECU 61 obtains a low temperature target EGR 
rate basic value MegrO by multiplying standard target 
EGR rate basic value MegrO and a correction coefficient 

15 0.2 as target EGR rate basis value MegrO (Megr0=Megr0 X 
0.2). When ATSstate=2, ECU 61 obtains NOx recovery 
target EGR rate basic value MegrO by multiplying standard 
target EGR rate basic value MegrO and a correction 
coefficient 0.8 as target EGR rate basis value MegrO 

20 (MegrO=MegrOXo.8) . When ATSstate = 3, ECU 61 obtains a 

desulf urization mode target EGR rate basic value MegrO by 
multiplying standard target EGR rate basic value MegrO 
and a correction coefficient 0 as target EGR rate basis 
value MegrO (Megr0=Megr0 X o ) . When ATSstate=4, ECU 61 

25 sets filter recovery mode target EGR rate basic value 

MegrO by multiplying standard target EGR rate basic value 
MegrO and a correction coefficient 0.5 as target EGR rate 
basis value MegrO (Megr0=Megr0 X o . 5 ) . 
[0129] When ATSstate=3 or 4 and when one of 

30 desulf urization mode and filter recovery mode is executed, 
if exhaust gas rising mode is selected, ECU 61 sets 
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exhaust gas rising mode target EGR rate basic value MegrO 
obtained by multiplying standard target EGR rate basic 
value MegrO and a correction coefficient 0 as target EGR 
rate basis value MegrO (MegrO=MegrO X o ) . Therefore, ECU 
5 61 stops EGR. When a malfunction avoiding mode is 

selected, ECU 61 sets malfunction avoiding mode target 
EGR rate basic value MegrO obtained by multiplying 
standard target EGR rate basic value MegrO and a 
correction coefficient 0,8 as target EGR rate basis value 
10 MegrO (MegrO=MegrO X o . 8 ) . 

[0130] Fig. 31 is a flowchart for a calculation 
routine of a target EGR rate, which is executed by module 
M5 of ECU 61. 

[0131] At step S551 ECU 61 reads target EGR rate basic 
15 value MegrO and intake system response time constant Kkin. 
At step S552 ECU 61 executes a delay processing of target 
EGR rate basic value MegrO using the following expression 
(19) which includes intake system response time constant 
Kkin, and stores the obtained value as Megrd. 

20 Megrd=Megrd n _i X ( l-Kkin) +Mger0 X Kkin (19) 

[0132] At step S553 ECU 61 calculates target EGR rate 
Megr by executing an advance processing of Megrd using 
the following expression (20) which employs GKeegr as a 
coefficient . 

25 Megr=Gkeegr X MegrO- (GKeegr- 1) X Megrd (20) 

[0133] Fig. 32 is a flowchart for a calculation 
routine of a target EGR gas quantity, which is executed 
by module M5 of ECU 61. 

[0134] At step S561 ECU 61 reads target intake air 
30 quantity tQac, target EGR rate Megr and intake system 

response time constant Kkin. At step S562 ECU 61 obtains 
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a target EGR gas quantity basic value QecO by multiplying 

trget intake air quantity tQac and target EGR rate Megr 

(QecO=tQacXMegr ) . At step S563 ECU 61 executes a delay 

processing of target intake air quantity tQac using the 
5 following expression (21) which includes intake system 

response time constant Kkin, and stores the obtained 

value as tQecd. 

tQecd=tQecd n -! X (l-Kkin) +tQec0 X Kkin (21) 

[0135] At step S564 ECU 61 obtains a target EGR gas 
10 quantity tQec by executing an advance processing of tQecd 

using the following expression (22) which includes intake 

system response time constant Kkin. 

tQec=GKqec X tQecO- (Gkqec-1 ) X tQecd (22) 

[0136] Fig. 33 is a flowchart of a calculation routine 
15 of a target EGR valve opening, which is executed by 

module M5 of ECU 61. 

[0137] At step S571 reads target EGR gas quantity tQec 
and EGR gas flow velocity Cqe . At step S572 ECU 61 
obtains a target EGR valve opening basic value tAegrO by 

20 diving target EGR gas quantity tQec by EGR gas flow 
velocity Cqe ( tAegrO=tQec/Cqe ) . At step S573 ECU 61 
calculates a target EGR valve opening tAegr from the 
following expression (23) on the basis of target EGR 
valve opening basis value t AegrO . The calculation of 

25 target EGR valve opening tAegr depends on a calculation 
method based on a Venturi model. 

tAegr=tAegrO/{ ^(1 -(t AegrO/ AEGRB# f) ( 23 ) 

where AEGRB# is a representative cross-sectional area of 
EGR passage. 
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[0138] At step S574 ECU 61 obtains an EGR valve step 
number STEPEGR by converting target EGR valve opening 
tAegr using a table shown in Fig. 34. 

[0139] Fig. 35 is a flowchart of a calculation routine 
5 of a target turbine nozzle opening, which is executed by 
module M5 of ECU 61. 

[0140] At step S581 ECU 61 reads engine speed Ne, 
target EGR rate Megr and target acceleration request 
injection quantity Qfdrv. At step S582 ECU 61 retrieves 

10 a turbine nozzle opening basic value TravO for achieving 
target excess coefficent Tlamb and target EGR rate Megr 
from a map shown in Fig. 36 with reference to engine 
speed Ne and target acceleration request injection 
quantity Qfdrv. At step S583 ECU 61 retrieves a turbine 

15 nozzle opening correction value Travq from a map shown in 
Fig. 37 with reference to engine speed Ne and target 
acceleration request injection quantity Qfdrv. At step 
S584 ECU 61 obtains a target turbine nozzle opening Trav 
by adding turbine nozzle opening basic value TravO and 

20 turbine nozzle opening correction value Travq 
(Trav=TravO+Travq) . 

[0141] Fig. 38 is a flowchart of a response delay 
compensation routine of a target turbine nozzle opening 
Trav, which is executed by module M5 of ECU 61. 

25 [0142] Variable nozzle type turbocharger 15 generates 
a response delay of gas flow and an operational delay of 
actuator 15 for driving a variable vane of turbine 15b. 
The response delay varies according to an exhaust gas 
flow rate Qexh on the assumption that the response delay 

30 includes operational delays of compressor 15a and turbine 
15b. The operational delay of actuator 51 is constant. 
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At steps S593 and S594 ECU 61 compensates the response 
delay, and at steps S596 and S597 ECU 61 compensates the 
operational delay. 

[0143] At step S591 ECU 61 reads target turbine nozzle 
5 opening Trav and exhaust gas flow rate Qexh. At step 
S592 ECU 61 retrieves an exhaust system response time 
constant Tcvgt from a table shown in Fig. 39 with 
reference to exhaust gas flow rate Qexh, and retrieves an 
advance compensation coefficient Gkvgt from a map (table) 

10 shown in Fig. 40 with reference to exhaust gas flow rate 
Qexh. At step S593 ECU 61 executes a delay processing of 
target EGR rate Megr using the following expression (24) 
which includes exhaust system response time constant 
Tcvgt and stores the obtained value as Travd. 

15 Travd=Travd n _i X ( 1-Tcvgt ) +Trav XTcvgt (24) 

[0144] At step S594 ECU 61 executes an advance 
processing of target turbine nozzle opening Trav using 
the following expression (25) which includes advance 
compensation coefficient GKvgt , and stores the obtained 

20 value as Travff . 

Travff=GKvgt X Trav- (GKvgt- 1 ) X Travd (25) 

[0145] At step S595 ECU 61 obtains a sum of Traveff 
and Travefb and stores the obtained value as Travc 
(Travc=Travf f +Travf b) , wherein Travfb is a feedback 

25 correction quantity obtained on the basis of target 

intake air quantity tQac and intake air quantity Qac. 
[0146] At step S596 ECU 61 executes a delay processing 
of Travc using the following expression (26) which 
includes a drive system response time constant TCACT# , 

30 and stores the obtained value as Travcd. 

Travcd=Travcd n _i X(i-TCACT#) XTravc XTCACT# (26) 
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[0147] At step S597 ECU 61 calculates a target turbine 
nozzle opening Travf by executing an advance processing 
of Travc using the following expression (27) which 
employs GKACT# as a coefficient. 

5 Travf =CKACT# X Travc- ( GKACT#- 1 ) XTravcd (27) 

[0148] Fig. 41 is a flowchart of a target duty ratio 
calculation routine which is executed by module M5 of ECU 
61. 

[0149] At step S601 ECU 61 reads target turbine nozzle 
10 opening Travf. At step S602 ECU 61 retrieves a target 

duty ratio VNduty, which is a signal of driving actuator 

51 from a map (table) shown in Fig. 42 with reference to 

target turbine nozzle opening Travf. 

[0150] Fig. 43 is a flowchart of a target main 
15 injection timing calculation routine, which is executed 

by module M5 of ECU 61. 

[0151] At steps S611 and S612 ECU 61 reads mode 
decision value ATSstate and retrieves a target main 
injection timing basic value MIT0 according to the 

20 operation mode from a map corresponding to mode decision 
value ATSstate. In this embodiment, when the recovery 
processing is executed, ECU 61 corrects target main 
injection timing basic value MIT0 retrieved from a normal 
map (standard) according to the target exhaust gas 

25 temperature, and sets the corrected value as a recovery 
mode target main injection timing basis value MIT0. 
Recovery mode target main injection basic value MIT0 is 
set at a timing retarded from a top dead center. 
[0152] More specifically, when ATSstate=l, ECU 61 

30 retrieves a normal mode target main injection timing MIT0 
from a reference map shown in Fig. 44. When ATSstate=2, 
ECU 61 sets NOx recovery mode target main injection 
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timing basic value MITO at a value obtained by retarding 
MITO of the reference mode by 10Q (crank angle) 
(MIT0=MIT0 + 10QCA) . When ATSstate=3, ECU 61 sets 
desulf urization mode target main injection timing basic 
5 value MITO at a value obtained by retarding MITO of the 
reference mode by 10Q (crank angle) (MIT0=MIT0+10QCA) . 
When ATSstate=4, ECU 61 sets filter recovery mode target 
main injection timing basic value MITO at a value 
obtained by retarding MITO of the reference mode by 10Q 

10 (crank angle) (MIT0=MIT0+10QCA) . 

[0153] When ATSstate=3 or 4 and when one of 
desulf urization mode and filter recovery mode is executed, 
if exhaust gas resing mode is selected, ECU 61 sets 
exhaust gas rising mode target main injection timing 

15 basic value MITO at a value obtained by retarding MITO of 
the reference mode by 10Q (crank angle) (MIT0=MIT0+10QCA) . 
If malfunction avoiding mode is selected, ECU 61 sets 
exhaust gas rising mode target main injection timing 
basic value MITO at a value obtained by retarding MITO of 

20 the reference mode by 62 (crank angle) (MIT0=MIT0+6QCA) . 
[0154] At step S613 ECU 61 reads intake system 
response time constant Kkin and obtains a target main 
injection timing MIT by executing a delay processing of 
MITO using the following expression (28) which includes 

25 intake system response time constant Kkin. 

MIT=MIT n _i X (1-Kkin)+MIT0 XKkin (28) 

[0155] At step S614 ECU 61 sets main injection timing 
MITf by adding target main injection timing MIT and a 
main injection timing correction value MITfb 

30 (MITf =MIT+MITfb) . When ATSstate=0, ECU 61 determines 

main injection timing MIT by executing a low temperature 
mode ignition timing control routine. 
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[0156] As discussed above, module M6 of ECU 61 rises 
the exhaust gas temperature by retarding the main 
injection timing and advances the pilot injection timing 
before the normal timing to suppress the generation of 
5 smoke and to reduce combustion noise. The pilot 

injection timing may be set in a manner as is similar to 
that of the main ignition timing. That is, a pilot 
injection timing basic value obtained from the normal 
mode map is advanced by a predetermined angle, and a 

10 delay processing of the obtained value is executed. 

[0157] There is discussed an air/fuel ratio feedback 
control. Fig. 4 5 shows a flowchart of an air/ fuel ratio 
feedback control routine, which is executed by module M6 
of ECU 61. Although the embodiment according to the 

is present invention has been shown and described to employ 
a PID algorithm represented by the following expression 
(29) of a proportion plus integral plus derivative 
compensator, the other algorithm may be employed. 

u(t) - Kp{e(t)+ — f e(t)dt + KD^^l + u(t0) - - - ( 29 ) 

\ KlJ dt J 

20 where u(t) is a manipulated variable, KP is a proportion 

gain, KI is an integral time constant, KD is a derivative 

time constant, e(t) is a difference, and u(t0) is an 

initial value. 

[0158] At step S701 ECU 61 reads target excess air 
25 ratio Tlamb, excess air ratio lamb, and mode decision 
value ATSstate. At step S702 ECU 61 calculates a 
disjunction (difference) dlamb between target excess air 
ratio Tlamb and excess air ratio lamb (dlamb = Tlamb - 
lamb) . 

30 [0159] At step S703 ECU 61 determines whether or not 
ATSstate=0, 2 or 3 . When the determination at step S703 
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is affirmative, the routine proceeds to step S704. When 
the determination at step S703 is negative, the routine 
proceeds to step S711. 

[0160] At each of steps S704 and S711, ECU 61 sets 
5 compensation gains KPlamb, Kllamb and KDlamb from tables 
shown in Figs. 46A, 46B and 46C, respectively, on the 
basis of excess air ratio lamb. At each of steps S705 
and S712, ECU 61 calculates an integral correction value 
Ilamb using the following expression (30). 

10 Ilamb=Ilamb n -i+(dT/KIlamb) X 6 lamb (30) 

[0161] At each of steps S706 and S713, ECU 61 limits a 
magnitude of integral correction value Ilamb within a 
predetermined range. At each of steps S707 and S714 # ECU 
61 calculates a derivative correction value Dlamb using 
15 the following expression (31). 

Dlamb=( 6 lamb- <5lamb n -i) XKDlamb/dT (31) 

[0162] At each of steps S708 and S715, ECU 61 
calculates a PID correction quantity Qffb, ETCfb (which 
includes a proportional term) from each of the following 
20 expressions (32A) and (32B). 

Qf f b=KPlmabf X ( 6 lamb+Ilambf +Dlambf ) +Klambf 0# 

(32A) 

ETCf bb=KPlmaba X ( 6 lamb+Ilamba+Dlamba) +Klamba0# 

(32B) 

25 where Klambf0# and Klambd0# are initial values of the 
respective correction values. 

[0163] At step S709 ECU 61 substitutes ETCfb n _i 
obtained in the previous routine in ETCfb (ETCfb = ETCfb n - 
i) . At step S710 ECU 61 calculates a final fuel injection 
30 quantity Qfdes by adding Qffb to target fuel injection 
quantity tQf (Qf des=Qf f b+tQf ) . 
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[0164] On the other hand, at step S716 ECU 61 
substitutes Qffb n _i obtained in the previous routine in 
Qffb (Qf fb = Qf fb n -i) . At step S717 ECU 61 calculates a 
final intake throttle value opening ETCf by adding ETCb 
5 to intake throttle valve opening ETC ( ETCf =ETC+ETCf b ) . 
[0165] There is discussed a calculation of a target 
exhaust gas temperature. Fig. 47 is a flowchart of a 
target exhaust gas temperature calculation routine, which 
is executed by module M8 of ECU 61. 

10 [0166] At steps S801 and S802 ECU 61 reads mode 

decision value ATSstate, selects a map corresponding to 
mode decision value ATSstate, and calculates a target 
exhaust gas temperature basic value tTexhO according to 
the selected map. That is, when ATSstate=3, ECU 61 sets 

15 a desulfurization mode target exhaust gas temperature 
basic value tTexhO at 730QC. When ATSstate=4, ECU 61 
sets a filter recovery mode target exhaust gas 
temperature basic value tTexhO at 6702C. Further, when 
an exhaust gas rising mode is selected, ECU 61 sets an 

20 exhaust gas rising mode target exhaust gas temperature 
basic value tTexhO at 700QC. 

[0167] At step S803 ECU 61 determines a target exhaust 
gas temperature tTexh by executing a delay processing of 
basic value tTexhO using the following expression (33) 

25 which includes intake system response time constant Kkin . 

tTexh=tTexh n _i X ( 1-Kkin ) + tTexhO X Kkin ( 33 ) 

[0168] There is discussed a calculation of the exhaust 
gas temperature. Fig. 48 is a flowchart of the exhaust 
gas temperature calculation routine, which is executed by 

30 module M8 of ECU 61. 

[0169] At step S901 ECU 61 calculates an output vTexh 
of sensor 73. At step S902 ECU 61 obtains exhaust gas 
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temperature TexhO by converting vTexh using a table shown 
in Fig. 49. At step S903 ECU 61 executes a weighted 
average processing of exhaust gas temperature TexhO and 
sets the obtained value as exhaust gas temperature Texh. 
5 [0170] There is discussed a feedback control of the 
exhaust gas temperature. Fig. 50 is a flowchart of the 
feedback control routine of the exhaust gas temperature, 
which is executed by module M9 of ECU 61. 
[0171] Although the embodiment according to the 
10 present invention has been shown and described to employ 
a PID algorithm represented by the following expression 
(34) of a proportion plus integral plus derivative 
compensator, the other algorithm may be employed. 

u(t) = Kp{e(t) + — fe(t)dt + KD^^l + u(t0) - - - ( 34 ) 

[ KlJ dt J 

15 where u(t) is a manipulated variable, KP is a proportion 

gain, KI is an integral time constant, KD is a derivative 

time constant, e(t) is a difference, and u(t0) is an 

initial value. 

[0172] At step S1001 ECU 61 reads target exhaust gas 
20 temperature tTexh and exhaust gas temperature Texh. At 
step S1002 ECU 61 calculates a disjunction (difference) 
dTexh between target exhaust gas temperature tTexh and 
exhaust gas temperature Texh (dTexh = tTexh - Texh) . At 
step S1003, ECU 61 determines proportion, integral and 
25 derivative compensation gains KPlamb, Kllamb and KDlamb 

from tables shown in Figs. 51A, 51B and 51C, respectively, 
on the basis of excess air ratio lamb. At step S1004, 
ECU 61 calculates an integral correction value I texh 
using the following expression (35). 
30 Itexh=Itexhn-l+(dT/KItexh) X <5 texh (35) 
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[0173] At step S1005, ECU 61 limits a magnitude of 
integral correction value Itexh within a predetermined 
range. At step S1006, ECU 61 calculates a derivative 
correction value Dtexh using the following expression 
5 (36). 

Dtexh=( 6 texh- 5 texh n . x ) XKDtexh/dT (36) 

[0174] At step S1006, ECU 61 calculates a PID 
correction quantity MITfb (which includes a proportional 
term) from the following expression (37). 

10 MITfb=KPtexhX ( 6 texh-Itexh+Dtexh) +Ktexh0# (37) 

where Ktexh0# is an initial value of the correction value. 
[0175] With the thus arranged embodiment according to 
the present invention, it becomes possible to derive the 
following advantages . 

15 [0176] During the desulf urization processing of NOx 
trap catalyst 32 and the filter recovery processing of 
diesel particulate filter 33, exhaust gas temperature 
Texh is risen to target temperature tTexh which is higher 
than the normal mode temperature, and air excess air 

20 ratio lamb is maintained at target excess air ratio tlamb 
according to the selected recovery mode. Therefore, even 
if the engine operating condition is changed due to the 
vehicle acceleration or if a traveling circumstance of 
the vehicle is changed, the system according to the 

25 present invention prevents excess air ratio lamb from 

changing according to these changes. This prevents the 
deterioration of NOx trap catalyst 32 and the generation 
of malfunction such that an element of diesel particulate 
filter 33 is cracked. 

30 [0177] Figs. 52A through 52C respectively show a 

relationship between the air/fuel ratio and the exhaust 
gas temperature, a relationship between the air /fuel 
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ratio and the CO discharge quantity, and a relationship 
between the air/fuel ratio and the HC discharge quantity. 
The CO discharge quantity and the HC discharge quantity 
are the quantity of carbon monoxide and the quantity of 
5 hydrocarbon which are discharged form engine 1 per unit 
time. During the desulf urization processing, the 
air/fuel ratio is set at the stoichiometric air/fuel 
ratio or rich state to decompose sulfur content trapped 
in NOx trap catalyst 32. Exhaust gas temperature has a 

10 characteristic that the exhaust gas temperature rises as 
the air/fuel ratio is decreased. Accordingly, when the 
air/fuel ratio becomes out of the target range due to the 
change of the engine operating condition, the exhaust gas 
temperature excessively rises, and therefore an excessive 

15 heat load may be applied to NOx trap catalyst 32. 
[0178] Further, when the air/fuel is set at a 
stoichiometric air/fuel ratio or rich state, the CO 
discharge quantity and the HC discharge quantity become 
increased. Therefore, under this control state, if the 

20 air/fuel ratio is largely increased to a value outside of 
a target range, a reduction agent such as carbon monoxide 
radically reacts in the catalyst, and therefore, an 
excessive heat load may be applied to NOx trap catalyst 
32. Generally, NOx trap catalyst 32 has a limitation in 

25 heat resistance, and it is difficult to improve this 
limitation. 

[0179] Fig. 53 shows a relationship between a heated 
temperature of NOx trap catalyst 32 and a NOx conversion 
ratio of the NOx trap catalyst 32 which has been put in 
30 the heated temperature. As is apparent from Fig. 53, if 

NOx trap catalyst 32 once receives an excessive heat load, 
the performance of the catalyst is largely deteriorated. 
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[0180] According to the present invention, during the 
desulf urization recovery mode, even if the engine 
operating condition is varied, excess air ratio lamb is 
maintained constant. Therefore, it becomes possible to 
5 prevent NOx trap catalyst 32 from receiving excessive 

heat load and thereby preventing the deterioration of the 
performance of NOx trap catalyst 32. Further, it is 
preferable that the desulf urization recovery mode target 
exhaust gas temperature is set at a value lower than or 
10 equal to 750QC, and the upper limit thereof is around 
800QC. 

[0181] Figs. 54A and 54B show relationships of the 
exhaust gas temperature and a particulate combustion 
speed relative to the air/fuel ratio. The particulate 

15 combustion speed is a decreased quantity per unit time of 
particulates deposited on diesel particulate filter 33. 
During the filter recovery processing, excess air ratio 
lamb is set at a lean state so as to suitably suppress 
the combustion of particulates. The particulate 

20 combustion speed largely varies according to the change 
of the air/fuel ratio and has a characteristic that the 
particulate combustion speed largely increases as the 
air/fuel ratio is increased. On the other hand, when the 
air/fuel ratio is decreased to a value outside of the 

25 target range due to the change of the engine operating 
condition, there is a possibility that excessive heat 
load is applied to diesel particulate filter 33 and 
therefore a filter element 331 generates a crack A as 
shown in Fig. 55B or loses stoppers 333 as shown by 

30 reference A in Fig. 55C. If the increased quantity of 

the fuel injection quantity is further large, there is a 
possibility that discharged fuel cools diesel particulate 
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filter 33 and prevents the recovery operation. However, 
according to the present invention, during the filter 
recovery mode, excess air ratio lamb is maintained 
constant, and this prevents diesel particulate filter 33 
from receiving excessive heat load and the recovery 
thereof from being prevented by such a cooling due to the 
excessive fuel increase . 

[0182] This application is based on Japanese Patent 
Application No. 2003-114717 filed on, April 18, 2003 in 
Japan. The entire contents of this Japanese Patent 
Application are incorporated herein by reference. 
[0183] Although the invention has been described above 
by reference to certain embodiments of the invention, the 
invention is not limited to the embodiments described 
above. Modifications and variations of the embodiments 
described above will occur to those skilled in the art, 
in light of the above teaching. The scope of the 
invention is defined with reference to the following 
claims . 



